Probiotic Supplementation for Promotion of Growth in Children: A Systematic Review and Meta-Analysis

Probiotics are commonly prescribed to promote a healthy gut microbiome in children. Our objective was to investigate the effects of probiotic supplementation on growth outcomes in children 0–59 months of age. We conducted a systematic review and meta-analysis which included randomized controlled trials (RCTs) that administered probiotics to children aged 0–59 months, with growth outcomes as a result. We completed a random-effects meta-analysis and calculated a pooled standardized mean difference (SMD) or relative risk (RR) and reported with a 95% confidence interval (CI). We included 79 RCTs, 54 from high-income countries (HIC), and 25 from low- and middle-income countries (LMIC). LMIC data showed that probiotics may have a small effect on weight (SMD: 0.26, 95% CI: 0.11–0.42, grade-certainty = low) and height (SMD 0.16, 95% CI: 0.06–0.25, grade-certainty = moderate). HIC data did not show any clinically meaningful effect on weight (SMD: 0.01, 95% CI: −0.04–0.05, grade-certainty = moderate), or height (SMD: −0.01, 95% CI: −0.06–0.04, grade-certainty = moderate). There was no evidence that probiotics affected the risk of adverse events. We conclude that in otherwise healthy children aged 0–59 months, probiotics may have a small but heterogenous effect on weight and height in LMIC but not in children from HIC.


Introduction
The role of gut microbiota in human health has been studied extensively in the recent past [1]. Observational and experimental studies from community settings have shown that gut microbiota immaturity or dysbiosis is associated with risk of development of acute malnutrition and linear growth failure in children [2][3][4][5]. Probiotics are one of the commonly used products to promote healthy gut microbiomes, encompassing a USD 54 billion industry around the globe [6,7]. Probiotics are defined as live microorganisms that, when administered in adequate amounts, confer a health benefit to the host [8]. Prebiotics are a non-digestible food ingredient that benefit the host by stimulating the growth or activity of microorganisms indigenous to the human digestive tract [9]. Synbiotics are a combination of both probiotics and prebiotics [10]. Multiple studies have assessed the usage of probiotics

Materials and Methods
We conducted a systematic review and meta-analysis and reported our findings according to Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. We developed a team consensus on the study questions and methods in a protocol before the start of the study and registered this protocol on PROSPERO (CRD42020154352).
We included individual-and cluster-randomized trials. Trials with multiple treatment arms were included if the only difference between the arms was probiotic usage. We included studies with children aged 0 to 59 months who were supplemented with oral probiotics or synbiotics. We excluded studies that only tested prebiotics. As our review considered both probiotics and synbiotics, the term "probiotic" hereafter refers to both probiotics and synbiotics. The studies were included irrespective of dose, dosage forms, and strain of probiotics. We included studies that had an adequate comparison group such as standard-of-care, placebo, or no-intervention groups. We excluded observational studies such as cohort studies, case-control studies, case series, and case reports. We excluded studies on children with congenital abnormalities, syndromic diagnosis, and chronic conditions such as cystic fibrosis or inflammatory bowel disease. Furthermore, we excluded studies in which the authors declared that participants were already undernourished or malnourished.
We searched multiple electronic databases, including PubMed, Embase, Cochrane Central Register of Controlled Trials, CINAHL, Scopus, and LILACS. The last date of search was 6 November 2020. The search strategy for all databases is available in the supplementary document (Text S1, supplementary document). We searched ClinicalTrials.gov for ongoing studies. We also searched the reference sections of published studies and systematic reviews. We contacted the authors in relevant fields for any new studies. If growth data were measured but not reported in the study paper, we contacted the authors for those results. Trial registries were searched and checked for unpublished data as well as newly published data. If a study was only available in a language other than English, we attempted to translate the paper. If an adequate translation was unavailable, the study was excluded.
At least two authors screened titles and abstracts using Covidence software [14], and extracted data in duplication (AI, JC, NP, JE, MZ ES, CF). We double extracted relevant data using a data collection form specifically designed for this review. Data were extracted independently into the same form. We extracted the data for study design, study setting (hospital vs. community, country, country income status), inclusion and exclusion criteria, participant characteristics (age, nutritional status, gender), and characteristics of intervention (type, strain, form, duration, frequency, dose, comparison group). We included data on all outcomes at the longest follow-up reported by individual studies. If longest follow-up was not reported, we used suitable data the study provided which was either in the form of total growth gain over the course of the study or growth gain per time. Any duplicate data were only counted once. Any disagreement among the authors during any stage of the study was resolved by discussion and review of the publication(s) with consultation of the senior author (AI).
Our primary outcomes included weight-for-age (continuous outcome, Z scores or kg) and height-for-age (continuous outcome, cm or Z scores). Growth data were used in the form of Z scores as to WHO standards or in the primary units. Secondary outcomes Nutrients 2022, 14, 83   3 of 18   included weight-for-height, BMI, head circumference, and adverse events such as nausea,  vomiting, diarrhea, abdominal pain, flatulence, and sepsis. We pooled the dichotomous outcomes to obtain a summary estimate in the form of mean relative risk (RR) and reported with its 95% confidence interval (CI). We used the standardized mean difference effect size for continuous outcomes due to studies reporting data in different units (e.g., a few studies reported weight in kg and the others in Z scores) and reported the standardized mean effect with its 95% CI. We used Review Manager 5.4 and Stata to conduct the meta-analysis [15,16]. We used the random-effects models to pool data as the effect of probiotics could be different in different study populations. We used funnel plots and Egger's test to assess for publication bias.
For eligible studies with multiple treatment arms, one eligible pair was selected and included, and if more than two groups were eligible, they were combined into a single pairwise comparison. If a trial had multiple arms that addressed different doses, those arms were combined and compared with the control arm to avoid double counting the control arm in the synthesis. Cluster trials were synthesized together with individually randomized trials using cluster adjusted values. If trial results were not cluster adjusted, we adjusted the result by methods given in the Cochrane handbook [17]. Extended details of data input are found in the supplementary document (Text S2).
We assessed clinical, methodological, and statistical heterogeneity of effects reported in the literature. Statistical heterogeneity was assessed using the χ 2 test, I 2 , and tau-squared statistics, and visual inspection of the forest plot. We considered statistical heterogeneity to be significant if the p value of the χ 2 was <0.1, I 2 values were above 50%, and forest plots showed different magnitude and effect of the intervention. We conducted subgroup analyses to explore reasons for any substantial statistical heterogeneity.
We conducted all analyses for low-and middle-income countries (LMIC) separately from high-income countries (HIC), given presumed heterogeneity in study populations across these settings. A subgroup analysis was conducted for growth outcomes using a χ 2 test to assess whether the effects of probiotics were significantly different for the following subgroups: age: 0-<6 months vs. 6-<60 months, probiotic interventions with single vs. multiple strains, synbiotics vs. probiotics, and participant status of healthy vs. premature/low birth weight. The healthy group included participants not defined as premature/low birth weight. As previously described, we excluded all studies in which authors declared that participants were undernourished or malnourished. We also completed a post hoc analysis in which we calculated the effect of each strain, or combination of strains, for each of our main outcomes.
We conducted sensitivity analyses by excluding studies with high risk of bias, those studies where data were supplied in a form other than mean (SD), or when effect-size data were extracted from figures.
Two authors assessed and agreed upon bias using the Cochrane Collaboration's riskof-bias tool-2 (ROB-2) for assessing risk of bias for all outcomes from the included studies in the meta-analysis [18]. Using this tool, results for each outcome were judged as either low, some concerns, or high risk of bias. The certainty of overall evidence for the effect of probiotics for an outcome was assessed using the Grading of Recommendations Assessment, Development and Evaluation (GRADE) method [19]. We present the results of the quality assessment in the form of a summary-of-findings table separately for high-income countries and low-and middle-income countries.

Literature Search
Our literature search identified 11,158 titles after exclusion of duplicates. Figure 1 shows the results of the literature search. The 11,158 studies were reduced to 243 full-text studies after applying the exclusion criteria stated in our methods section. After screening the full text of 243 studies, we ultimately included 79 studies in our systematic review, of which 54 studies were from high-income countries, and 25 studies from low-and middle- income countries . We excluded 164 studies, and reasons for exclusion can be found in Table S1. , x FOR PEER REVIEW 4 of 18 the full text of 243 studies, we ultimately included 79 studies in our systematic review, of which 54 studies were from high-income countries, and 25 studies from low-and middleincome countries . We excluded 164 studies, and reasons for exclusion can be found in Table S1.

Weight-for-Age
In LMIC, twenty-one studies reported data on weight-for-age and included a total of 8417 participants (4323 probiotics, 4094 control) [21,26,27,33,36,38,39,50,58,60,61,[68][69][70]72,77,[83][84][85]91,92]. The results showed low-certainty evidence that probiotics had a small effect on weight when compared to the control group (SMD: 0.26, 95% CI: 0.11-0.42, p = 0.001, I 2 = 87%, Figure 2). The GRADE evidence was downgraded to low due to statistical heterogeneity in the pooled data and clinical heterogeneity in the use of probiotics in the included studies ( Table 1). The funnel plot and Egger's test showed no evidence of publication bias (p = 0.07, Funnel plot in Figure S1). The summary risk of bias for weightfor-age in low-and middle-income countries is shown in Figure 3. The results were similar without high-risk-of-bias studies (SMD 0.31, 95% CI: 0.13-0.48, p < 0.001, I 2 = 89%). A subgroup analysis comparing single versus multiple strain interventions found that singlestrain interventions had a greater impact on weight gain than multiple strain interventions (P subgroup = 0.002). Figure 4 shows the results of other sensitivity analyses and subgroup analyses. All other subgroup and sensitivity analyses did not show significantly different results for weight-for-age. A post hoc subgroup analysis based on type of strains did not show a particular individual strain or combination of strains that were effective in terms of effect on weight (Table S6). effect of probiotics varied in magnitude. We also downgraded for heterogeneity in the use of probiotics used in the included studies. c. All the studies were conducted in low-r and middle-income countries d. Overall sample size from all the included studies in the meta-analysis was more than 8000. The CI did not include 0. e. None of the included studies in this analysis were at high risk of bias. f. Even though the statistical heterogeneity was only 25%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. g. The I 2 was 57% h. The overall sample size was less than 2000 and the confidence interval of the summary estimate included 0. i. Even though the statistical heterogeneity was only 0%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. j. The overall sample size of the included studies was less than 400 and the confidence interval of the summary estimate was wide and included 0. k. Even though one of the included studies was at high risk of bias for this outcome, a sensitivity analysis by excluding this study did not change the magnitude, direction, or statistical significance of the summary estimate. l. Even though the statistical heterogeneity was only 20%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. Abbreviations: CI: confidence interval; SMD: standardized mean difference; RR: risk ratio, BMI: body mass index RCT: randomized controlled trial.

MODER-ATE
Footnotes: a. Even though three of the included studies were at high risk of bias for this outcome, a sensitivity analysis by excluding these studies did not change the magnitude, direction, or statistical significance of the summary estimate. b. The I 2 was 87%. Inspection of the forest plot showed the effect of probiotics varied in magnitude. We also downgraded for heterogeneity in the use of probiotics used in the included studies. c. All the studies were conducted in low-r and middle-income countries d. Overall sample size from all the included studies in the meta-analysis was more than 8000. The CI did not include 0. e. None of the included studies in this analysis were at high risk of bias. f. Even though the statistical heterogeneity was only 25%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. g. The I 2 was 57% h. The overall sample size was less than 2000 and the confidence interval of the summary estimate included 0. i. Even though the statistical heterogeneity was only 0%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. j. The overall sample size of the included studies was less than 400 and the confidence interval of the summary estimate was wide and included 0. k. Even though one of the included studies was at high risk of bias for this outcome, a sensitivity analysis by excluding this study did not change the magnitude, direction, or statistical significance of the summary estimate. l. Even though the statistical heterogeneity was only 20%, we downgraded for clinical heterogeneity in the use of probiotics in the included studies. Abbreviations: CI: confidence interval; SMD: standardized mean difference; RR: risk ratio, BMI: body mass index RCT: randomized controlled trial.

Height-for-Age
Twelve studies from LMIC reported data on height-for-age and included a total of 2561 participants (1347 probiotics, 1214 control) [21,27,33,36,39,50,70,[83][84][85]91,92]. The results provided moderate-certainty evidence that probiotics had a small effect on heightfor-age when compared to the control group (SMD: 0.16, 95% CI: 0.06-0.25, p = 0.002, I 2 = 25%, Figure 5). The GRADE evidence was downgraded due to heterogeneity in the clinical use of probiotics in the included studies ( Table 1). The funnel plot and Egger's test showed no evidence of publication bias (p = 0.25). All subgroup and sensitivity analyses did not show significantly different results for height-for-age (Figure 4). A summary risk of bias for height-for-age and other outcomes in LMIC are available in Figures S6-S9 of the supplementary document. 25%, Figure 5). The GRADE evidence was downgraded due to heterogeneity in the clinical use of probiotics in the included studies ( Table 1). The funnel plot and Egger's test showed no evidence of publication bias (p = 0.25). All subgroup and sensitivity analyses did not show significantly different results for height-for-age (Figure 4). A summary risk of bias for height-for-age and other outcomes in LMIC are available in Figures S6-S9 of the supplementary document.

Other Outcomes
Among studies from LMIC, probiotics did not have any significant effect on other outcomes including adverse events (forest plots in Figures S10-S15). A funnel plot for head circumference in LMIC is available in Figure S3. The other analyses for LMIC, including subgroup and sensitivity analyses, are displayed in Figure 4.

Weight-for-Age
Fifty-one studies from HIC reported data on weight-for-age and included a total of 10,832 participants (5759 probiotics, 5074 control) [ . The results provided moderate-certainty evidence that probiotics did not have a clinically meaningful effect on weight when compared to the control group (SMD: 0.01, 95% CI: −0.04-0.05, p = 0.78, I 2 = 7%, Figure  S16). The GRADE evidence was downgraded due to clinical heterogeneity in the use of probiotics in the included studies ( Table 2). The funnel plot and Egger's test showed no evidence of publication bias (p = 0.64, Figure S2). The sensitivity analysis without highrisk-of-bias studies showed similar results (SMD 0.00, 95% CI: −0.04-0.05, p = 0.88, I 2 = 6%). A post hoc subgroup analysis based on type of strain did not show a strain or combination of strains that influenced any of the growth outcomes from studies from high-income countries (Table S7). All subgroup and sensitivity analyses did not show significantly different results for weight-for-age ( Figure S32).  Figure 5. Effect of use of probiotics on height-for-age in children 0-59 months of age from low-and middle-income countries [21,27,33,36,39,50,70,[83][84][85]91,92].

Other Outcomes
Among studies from LMIC, probiotics did not have any significant effect on other outcomes including adverse events (forest plots in Figures S10-S15). A funnel plot for head circumference in LMIC is available in Figure S3. The other analyses for LMIC, including subgroup and sensitivity analyses, are displayed in Figure 4.

Height-for-Age
Thirty-two studies from HIC reported data on height-for-age and included a total of 6118 participants(3350 probiotics, 2768 control) [20,22,25,28,30,31,37,40,41,44,45,49,[53][54][55][56]63,66,67,71,73,75,[78][79][80][86][87][88]90,94,96,97]. The results provided moderate evidence that probiotics did not have a clinically meaningful effect on height-for-age when compared to the control group (SMD: −0.01, 95% CI: −0.06-0.04, p = 0.71, I 2 = 0%, Figure S17). The GRADE evidence was downgraded due to clinical heterogeneity of probiotics used in the included studies. The funnel plot and Egger's test showed no evidence of publication bias (p = 0.87, Figure S4). All subgroup and sensitivity analyses did not show significantly different results for height-for-age ( Figure S32). Explanations: a. Even though four of the included studies were at high risk of bias for this outcome, the sensitivity analysis by the exclusion of these studies did not change the magnitude, direction, or statistical significance of the analysis. b. Even though the statistical heterogeneity was very low with I 2 values of 7%, we downgraded for clinical heterogeneity in the type of probiotics used in the included studies. c. All the studies were conducted in high-income countries. d. The confidence interval of the effect size includes 0, and we think the overall effect is about 0. The confidence intervals are narrow enough that we do not think that the summary size is imprecise. e. Even though the statistical heterogeneity was homogenous with I 2 values of 0%, we downgraded for clinical heterogeneity in the type of probiotics used in the included studies. f. Even though three of the included studies were at high risk of bias for this outcome, the sensitivity analysis by the exclusion of these studies did not change the magnitude, direction, or statistical significance of the analysis. g. The I 2 was 85% h. The total sample size of all the studies included in the meta-analysis was less than 1000. The CIs were wide i. Even though one of the included studies was at high risk of bias for this outcome, the sensitivity analysis by the exclusion of this study did not change the magnitude, direction, or statistical significance of the analysis. j. Even though the statistical heterogeneity was low with I 2 values of 26%, we downgraded for clinical heterogeneity in the type of probiotics used in the included studies. k. The total sample size of all the studies included in the meta-analysis was less than 1000. The CIs were wide, and increased risk cannot be excluded. Abbreviations: CI: confidence interval; SMD: standardized mean difference; RR: risk ratio, BMI: body mass index RCT: randomized controlled trial.

Other Outcomes
Probiotics did not show any significant effect on other primary and secondary outcomes in studies from HIC (forest plots in Figures S16-S25). The summary risk of bias for effect of probiotics for weight-for-age and other outcomes is shown in Figures S26-S31 of the supplementary document. A funnel plot for head circumference in HIC is found in Figure S5. The other analyses for HIC, including subgroup and sensitivity analyses, are displayed in Figure S32.

Discussion
This comprehensive systematic review evaluated the effects of probiotics on growth in children 0 to 59 months of age. Overall, there was no evidence that probiotics had a clinically meaningful effect on any of the growth outcomes in children from high-income countries. The data from low-and middle-income countries showed that there may be a small beneficial effect on weight and height gain; however, the certainty of evidence was low and moderate for these outcomes. There was no evidence that probiotics increased the risk of any of the adverse events including risk of sepsis from HIC and LMIC.
We used the GRADE approach to assess the overall certainty of evidence for the effect of probiotics on primary outcomes and selected secondary outcomes. The GRADE method of certainty assessment gives ratings of evidence for each outcome and considers factors such as type of study, risk of bias, inconsistency of results, indirectness of evidence, imprecision of the summary estimate, and publication bias [19]. All included studies were randomized trials, and we did not downgrade the certainty rating for study designs for any of the outcomes graded. Overall, there were a few studies that were at high risk of bias, and sensitivity analyses removing these studies did not substantively change the results, so we did not adjust the overall certainty of evidence rating for risk of bias. However, we did adjust the certainty rating due to clinical heterogeneity in the use of probiotics for all the outcomes graded. We conducted separate GRADE assessments for studies from highincome countries and those from low-and middle-income countries because we believe that environmental factors, diet, and gut microbiome might be different in these settings [2]. We therefore did not downgrade the certainty of evidence for indirectness for any of the graded outcomes. The number of studies that contributed data for an outcome varied among the outcomes. We downgraded the evidence for imprecision where the number of included studies was small and the confidence interval included a null effect.
Even though there was clinical heterogeneity in the use of probiotics, the pooled results were mostly homogenous around the null effect, especially from studies from highincome countries. Therefore, it can be concluded with reasonable confidence that probiotic supplementation in otherwise healthy children from high-income countries does not give any differential effect in terms of growth. The findings from trials conducted in low-and middle-income countries were mixed. The pooled results for the effects of probiotics on weight-for-age from low-and middle-income countries showed a small effect in favor of the intervention (SMD 0.26, 95% CI 0.11-0.42); however, there was significant statistical heterogeneity in the pooled data (I 2 = 87%) and the positive effect can be explained by three studies [27,36,92] that had an effect SMD > 1. Removal of these studies with outlying findings yielded results that were similar to those from high-income countries (SMD 0.05, 95% CI −0.02, 0.12). We adjusted the GRADE evidence by lowering the certainty grade to 'low' for this outcome, which means that we have low confidence in this estimate and future research might change this estimate.
Were there any subgroups that might be affected differently? The data from highincome countries did not show any differential effects of probiotics for groups such as age < 6 months vs. 6-59 months, probiotics vs. synbiotics, or single vs. multiple strains and nutritional status ( Figure S32). Thus, there does not seem to be an overall effect or any significant effect in subgroups from studies from high-income countries. The subgroup analyses from low-and middle-income countries are hard to interpret as the number of studies in each subgroup varied and some of the differences observed could be due to the small number of included studies. Nonetheless, there was evidence that single-strain probiotics may yield a more pronounced effect on weight-for-age compared to multiplestrain probiotics (Figure 4), although no particular single strain can be attributed to this result. This result was only present when pooling many studies that utilized a single-strain intervention. Future large studies will be required to further elaborate on and potentially replicate this finding. Indeed, more-targeted therapies could be more beneficial rather than using a single strain or combination of strains of probiotics [3].
This study is one of the largest systematic reviews conducted on the subject. We searched multiple databases and examined 11,158 titles and abstracts, and this included both published and ongoing studies. We did not apply any limitations of the literature search and did not exclude studies at the title/abstract screenings stages if they did not report outcomes in the abstract. This might be the reason that we were able to include many more studies compared to the last review published in 2014 [12]. We specified our analyses a priori and registered our protocol on a publicly available website before the review started. The major post hoc decision was to include studies from high-income countries as we understand that probiotics are commonly used to promote healthy gut microbiome in high-income countries, and it is important to review the available literature for their effect in promotion of healthy growth in children. We therefore conducted separate meta-analyses and certainty ratings for high-income countries and low-and middle-income countries. The limitations of the data presented in this study were that there was significant clinical heterogeneity in the use of probiotics with respect to type, duration, and combination. It is debatable if a meta-analysis should have been performed in the presence of such clinical heterogeneity. We believe that a meta-analysis was appropriate here as it served the purpose of assessing the overall magnitude and direction of effect from included studies. We adjusted the GRADE ratings for each outcome for clinical heterogeneity and conducted a post hoc subgroup analysis examining each strain or combination of strains and found no differential effect for any single strain or combination of strains for any of the growth outcomes assessed. Another potential limitation is related to type of growth outcomes included. Although we included a range of outcomes for growth, most of the outcomes were continuous and we did not include dichotomous outcomes such as undernutrition and stunting. We posit that if probiotics had any significant effect on prevention of stunting and undernutrition then that pattern would be mirrored in their effects on average height and weight gain.
There is a growing appreciation for the role of the gut microbiota, and accumulating evidence suggests an association between immaturity of the microbial community and undernutrition [4,99,100]. Although our understanding of the dynamic interplay between undernutrition and the microbiota is improving, empirical work examining therapeutic interventions to ameliorate microbiota dysbiosis remains limited. Methods that harness the microbiota through rationally designed microbial consortia, or manipulate the community through microbiota-directed therapeutic foods, show promise for treatment [3,4]. Our observation that children in low-and middle-income countries might benefit from probiotic supplementation is consistent with recent findings regarding the beneficial effects of a microbiota-targeted food intervention in Bangladesh [3]. There is a need for large-scale clinical trials that address the multifaceted role of the microbiota in childhood nutrition.

Conclusions
Probiotic supplementation does not seem to have a clinically meaningful effect on growth for apparently healthy children in high-income countries. However, there might be a small effect on weight and height in apparently healthy children from low-and middleincome countries. Future large-scale clinical trials are needed that assess the targeted therapies to prevent the gut dysbiosis associated with childhood undernutrition in lowand middle-income countries.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/nu14010083/s1, Text S1: Search strategies for databases through November 6, 2020; Text S2: Extended data-input details; Text S3: Results: Study Locations; Table S1: Excluded studies and reasons for exclusion; Table S2: Participant characteristics in the included studies from low-and middleincome countries; Table S3: Participant characteristics in the included studies from high-income countries; Table S4: Treatment and comparison characteristics in the included studies from lowand middle-income countries; Table S5: Treatment and comparison characteristics in the included studies from high-income countries; Table S6: Pooled analysis based on probiotic strains: Low-and middle-income countries; Table S7: Pooled analysis based on probiotic strains: High-income countries; Figure S1: Funnel Plot: Effect of probiotics on growth: Low-and middle-income countries: Weight; Figure S2: Funnel Plot: Effect of probiotics on growth: High-income countries: Weight; Figure S3: Funnel Plot: Effect of probiotics on growth: Low-and middle-income countries: Head Circumference; Figure S4: Funnel Plot: Effect of probiotics on growth: High-income countries: Height; Figure S5: Funnel Plot: Effect of probiotics on growth: High-income countries: Head circumference; Figure S6: Risk of bias in the included studies from low-and middle-income countries for effect of probiotics on the outcome of weight-for-age; Figure S7: Risk of bias in the included studies from low-and middle-income countries for the outcome of height-for-age; Figure S8: Risk of bias in the included studies from low-and middle-income countries for the outcome of head circumference; Figure S9: Risk of bias in the included studies from low-and middle-income countries for the outcome of sepsis; Figure S10: Forest Plot: Effect of use of probiotics on head circumference in children 0-59 months of age from low-and middle-income countries; Figure S11: Forest Plot: Effect of use of probiotics on BMI in children 0-59 months of age from low-and middle-income countries; Figure S12: Forest Plot: Effect of use of probiotics on weight for height in children 0-59 months of age from low-and middle-income countries; Figure S13: Forest Plot: Effect of use of probiotics on adverse events: Sepsis, in children 0-59 months of age from low-and middle-income countries; Figure S14: Forest Plot: Effect of use of probiotics on adverse events: Vomiting, in children 0-59 months of age from low-and middle-income countries; Figure S15: Forest Plot: Effect of use of probiotics on adverse events: Diarrhea, in children 0-59 months of age from low-and middle-income countries; Figure S16: Forest Plot Effect of use of probiotics on weight-for-age in children 0-59 months of age from high-income countries; Figure S17: Forest Plot: Effect of use of probiotics on height-for-age in children 0-59 months of age from highincome countries; Figure S18: Forest Plot: Effect of use of probiotics on weight-for-length in children 0-59 months of age from high-income countries; Figure S19: Forest Plot: Effect of use of probiotics on head circumference for age in children 0-59 months of age from high-income countries; Figure S20: Forest Plot: Effect of use of probiotics on BMI in children 0-59 months of age from high-income countries; Figure S21: Forest Plot: Effect of use of probiotics on adverse events: Sepsis, in children 0-59 months of age from high-income countries; Figure S22: Forest Plot: Effect of use of probiotics on adverse events: Vomiting, in children 0-59 months of age from high-income countries; Figure S23: Forest Plot: Effect of use of probiotics on adverse events: Diarrhea, in children 0-59 months of age from high-income countries; Figure S24: Forest Plot: Effect of use of probiotics on adverse events: Abdominal pain, in children 0-59 months of age from high-income countries; Figure S25: Forest Plot: Effect of use of probiotics on adverse events: Flatulence, in children 0-59 months of age from high-income countries; Figure S26: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of weight-for-age; Figure S27: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of height-for-age; Figure S28: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of weight-for-length; Figure S29: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of head circumference; Figure S30: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of BMI; Figure S31: Risk of bias in the included studies from high-income countries for effect of probiotics on the outcome of sepsis; Figure S32: Subgroup and Sensitivity analyses in high-income countries.